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ABSTRACT 

Analysis for the turbulent flow field and the 
effect of freestream turbulence on the surface heat 
transfer rate of a stagnation flow is presented. The 
emphasis is on the modeling of turbulence and its aug- 
mentation of surface heat transfer rate. The flow 
field considered is the region near the forward stag- 
nation point of a circular cylinder in a uniform 
turbulent mean flow. The freestream is steady and 
incompressible with a Reynolds number of the order of 
10 b and turbulence intensity of <5 percent. The flow 
field is divided into three regions: (1) a uniform 

freestream region where the turbulence is homogeneous 
and isotropic; (2) an external inviscid flow region 
where the turbulence is distorted by the variation of 
the mean velocity; and (3) an anisotropic turbulent 
boundary layer flow region. The turbulence modeling 
techniques are the k-e two-equation model in the 
external flow and the time-averaged turbulence trans- 
port equation in the boundary layer. The turbulence 
double correlations, the mean velocity, and the mean 
temperature within the boundary layer are solved 
numerically from the transport equations. The surface 
heat transfer rate is calculated as functions of the 
freestream turbulence longitudinal microlength scale, 
turbulence intensity, and Reynolds number. 


NOMENCLATURE 

coefficients of Eqs. (Al) and (A2), 
i = 1 to 5 mean velocity parameter 

B coefficients of Eq. (A5), i = 1 to 4 

c-j empirical constants, i = 1 to 5 

c p empirical constant 

D cylinder diameter 

E dimensionless turbulence dissipation rate 


F dimensionless mean flow streamwise velocity 

fi parametrical functions of <p 0 , i = 1 to 3 

G flow variable of Eqs. (Al) and (A2) 

H dimensionless mean temperature 

K dimensionless turbulence kinetic energy 

k turbulence kinetic energy 

i spatial distance 

Nu Nusselt number 

Pr Prandtl number, 0.7 

R sum of dimensionless Reynolds normal 

stresses, R x + Ry + R z 

Re Reynolds number 

R x dimensionless Reynolds normal stress along 

x-di recti on 

Ry dimensionless Reynolds normal stress along 

y-direction 

R z dimensionless Reynolds normal stress along 

z-direction 

S dimensionless Reynolds shear stress 

T mean temperature 

Tu freestream turbulence intensity 

t temperature fluctuation 

U,V mean flow velocity components 

v't transformed velocity of V, Eq. (11) 

u, v,w velocity fluctuations 
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x,y,z physical coordinate system 

6 stagnation point laminar boundary layer 

thickness 

e turbulence dissipation rate 

n transformed coordinate along y-direction 

e dimensionless turbulence correlation, 

a turbulence modeling length scale 

x turbulence longitudinal microlength scale 

v kinematic viscosity 

K transformed coordinate along x-direction 

p density 

<p azimuth angle 

time-averaged quantity 
Subscripts : 

0 condition near stagnation point (cp = 0.04°) 

1 condition at location 1 

2 condition at location 2 of boundary layer 
edge 

3 condition at location 3 of boundary layer 
edge 

e external or boundary layer edge condition 

0,e boundary layer edge condition at initial 

station 

F freestream condition 

m grid point along r or n direction 

w surface condition 


INTRODUCTION 

The momentum and thermal flow fields near the 
forward stagnation point of a circular cylinder in tur- 
bulent flow have been the focus of considerable study. 
When the turbulence intensity was 3 to 5 percent, 
existing experiments have shown that the local Nusselt 
number can increase by about 75 percent as compared 
with laminar flow conditions. With the current empha- 
sis on high turbine inlet temperature, it has become 
even more important to understand the effect of turbu- 
lence on the surface heat transfer, particularly in 
the leading edge stagnation region. To estimate this 
surface heat transfer rate, correlations derived from 
experimental data may be used. Detailed experimental 
data can be found in the paper by Lowery and Vachon 
U). An alternative is to rely on theoretical analysis 
to predict the surface heat transfer rate. The objec- 
tive of the present study is to develope an analysis 
which predicts the effect of turbulence on surface 
heat transfer rate. 

Existing boundary layer flow analyses (_2,3) indi- 
cated a large effect of turbulence on the stagnation 
point heat transfer rate. Other studies (4 to 8) 


showed that turbulence amplification because of vortex 
stretching can also influence the thin boundary layer, 
and therefore the surface heat transfer rate. Hijikata 
et al. (9J performed a theoretical and experimental 
study of the stagnation point anisotropic turbulence. 
This analysis showed that the turbulence length scale 
also affects the surface heat transfer rate. Traci 
and Wilcox (10j_ analytical ly studied the turbulence 
effect on the stagnation point flow field. The mean 
velocity gradients within the inviscid flow were found 
to be important in analyzing the stagnation point heat 
transfer. A difficulty in this type of analysis is 
obtaining a proper match of the boundary conditions 
between the different flow regions. Strahle (11 ) 
tried to resolve the turbulence matching conditions. 
Analytical relations between the turbulence properties 
at the stagnation point and in the freestream were 
found. In a previous study (12J Wang proposed a theo- 
retical boundary layer flow analysis and numerical 
computational procedure to investigate the turbulence, 
momentum, and thermal fields within the stagnation flow 
region. The numerical calculation showed that rapid 
amplification of the Reynolds normal stress along the 
boundary layer edge was required to predict approxi- 
mately the surface heat transfer rate as reported in 
existing experiments. 

Based on the experimental results and analytical 
methods of the existing studies (1^, to 12) , it is 
possible to pursue a detailed theoretical and numeri- 
cal analysis of the forward stagnation flow field of 
a circular cylinder in a turbulent freestream. 

Reported herein is such an analysis which predicts 
the surface heat transfer rate as functions of the 
freestream turbulence longitudinal microlength scale, 
turbulence intensity, and Reynolds number. The empha- 
sis is on the modeling of turbulence and its augmen- 
tation of the surface heat transfer rate. 

The flow field is divided into three regions: 

(1) a region away from the cylinder where the mean 
velocity is uniform and the turbulence is assumed to 
be homogeneous and isotropic, (2) an external inviscid 
flow region which has a constant rate of mean velocity 
variation and the mean velocity distorts the turbu- 
lence, and (3) a boundary layer region over the cylin- 
der surface. For turbulence modeling, the k-e 
two-equation modeling and the time-averaged turbulence 
transport equation are used respectively in Flow 
Regions 2 and 3. These turbulence modeling equations 
and the mean flow conservation equations are used to 
formulate the theoretical analysis. The turbulence 
double correlations, the mean velocity, and the mean 
temperature within the boundary layer region are 
solved numerically from the turbulence transport equa- 
tions and mean flow conservation equations. The solu- 
tion is obtained by matching the turbulence kinetic 
energy along the boundary between the flow regions. 
These analytical results are compared with the mea- 
surements from the existing experiments. 

ANALYTICAL FORMULATION OF FLOW REGIONS 

The flow region of interest is around the forward 
stagnation point of a circular cylinder in a turbulent 
freestream. Figure 1 schematically represents the flow 
field. The diameter of the cylinder is D and its axis 
is perpendicular to the freestream mean velocity Vp. 

The freestream Reynolds number is defined as 
Re = VpD/v, where \> is the kinematic viscosity. The 
velocity fluctuation in the freestream is vf and 
the freestream turbulence intensity is defined as 
Tu =y?p/V p. The freestream turbulence longitudinal 
microlength scale is denoted by xp . The cylinder 
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is assumed to have a constant surface temperature 
T w . The freestream has a mean temperature Tf. 

The stagnation point is taken as the coordinate- 
system origin with x and y coordinates parallel to and 
normal to the cylinder surface and z-coordinate par- 
allel to the cylinder axis. The mean flow field is 
assumed steady, incompressible, and two-dimensional. 
The mean velocity components are U and V along the 
x- and y coordinates. The turbulence is three- 
dimensional. The velocity fluctuations are denoted by 
u, v, and w along the x-, y-, and z-directions. The 
governing equations of the mean flow and turbulence 
are derived for each flow region. Following the 
laminar flow analysis, the relations 


c( 


f U e 

x)=| ^dx. 


U e y 


and n(x,y) = 


(2K) 


UTS 


( 1 ) 


are used for coordinate transformation. The analytical 
formulation, together with the assumptions and the 
boundary conditions, is described here. 


Region 1 

This flow region is far away from the cylinder. 
The flow is turbulent and its mean velocity is uni- 
form. The properties of the mean flow and its turbu- 
lence are not influenced by the existence of the 
cylinder. Appropriate conditions of the turbulence 
kinetic energy k, and the turbulence dissipation 
rate e, are required to characterize the turbulence. 
The turbulence is assumed to be homogeneous and iso- 
tropic. The freestream turbulence kinetic energy and 
the turbulence dissipation rate are given in Ref. 13 
as 
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3Vp 30v V p 

kp = 2 ‘ » and €p = 2 — (2) 

X F 

The turbulence longitudinal microlength scale is 
determined from the two point turbulence correlations 
(13) . The longitudinal two point turbulence correla- 
tion g, the spacial distance £, and the turbulence 
longitudinal microlength scale are related by the fol- 
lowing relations for small £. 


gU) 


(3) 


Region 2 

In this region, the mean velocity is disturbed by 
the cylinder. The mean velocity variation distorts 
the turbulence. It is assumed that the turbulence 
does not influence the mean velocity. Thus, the mean 
velocity is represented with the two-dimensional 
inviscid flow solution near the stagnation point. The 
mean velocity components are 

U = ax and V = -ay (4) 

with a = 4Vp/D. The turbulence kinetic energy 
variation along the stagnation point streamline is of 
interest in this study. The turbulence is modeled 
with the k-e two-equation turbulence modeling (14) . 


Introducing Eqs. (1) and (4) into the standard 
k and e equations, neglecting the diffusion terms 
in the x-direction and seeking the following forms of 
solutions: 


k = vaK(n) and g = va^E(v) 
the k and e equations become 

dK . K 2 . c d / K 2 dK 
n dn C ii E dn l C u c^E dn 

. dE _ _ 4c K + r 1 - L !<1 dE 

n d n " 4 1 y 2 k dn l C y c 4 E dn 


(5) 

( 6 ) 

(7) 


along the stagnation point streamline. The empirical 
constants are cl = 1.44, c 2 = 1.92 c 3 = 1.00, 

C 4 = 1.30, and c„ = 0.09. 

Equations ( 6 ) and (7) are solved for the turbulence 
properties K and E in the interval nl > n > n 2 » 
where nl is the boundary between Flow RegTons 1 
and 2, and n 2 is the boundary between Flow 
Regions 2 and 3. 

The boundary conditions are defined in the 
following: 

1. At n - nl, 

kp ep 

K = K i , and E = E 1= ^ ( 8 ) 


where kp and cp are previously determined by 
Eq. ( 2 ). 


2. Neglecting the diffusion term in Eq. (7) at 
nl, this equation gives the following additional 
boundary condition at n = n] . 


^ = - Lp - 4 C |C K, 

d n n 1 \ 2 K 1 1 u 1 j 

3. It is also assumed (1(D) that 
dK n 

d? = 0 at n = n 2 


(9) 


( 10 ) 


The finite difference method is used to solve the 
above boundary value problem. A summary of the com- 
putational procedure is given in Appendix A. The 
relation nl = Vp/(a v ) 0,5 is found from Eqs. (1) and 
(4), and n2 is resolved in the process of the 
numerical calculation. The computational results of 
the turbulence properties K and E at n = n2 are 
used to analyze the boundary layer flow region. 

Region 3 

A steady, two-dimensional, and incompressible 
turbulent boundary layer is assumed in the boundary 
layer region. The mean flow properties are described 
by the boundary layer continuity, momentum, and 
enthalpy equations. The momentum and enthalpy equa- 
tions contain turbulence double correlations, Uv and 
vt. A method to model the turbulence double correla- 
tions within a turbulent flow was developed by 
Donaldson et al. U5 and 16 ) . This method is used to 
formulate a theoretical analysis of the turbulence in 
this flow region. The time-averaged turbulence trans- 
port equations are derived following the theory in 
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Ref. 15. Assumptions are made for the turbulence 
closure relations and a modeling length scale A. 
Details of the theoretical formulation and assumptions 
can be found in Refs. 15 to 17. 

The mean flow velocity U e = ax is assumed along 
the boundary layer edge between regions 2 and 3. The 
mean flow properties and turbulence double correla- 
tions are also nondimensional ized as 


F -iL 

U e 


H = 


T T “2 

D _ JL_ R _ V_ D - W_ 

K v “ 9 » ~ 9 


u; 


u; 


u 


U e 


vt 


e e 


Introducing the above dimensionless variables and 
Eq. (1) into the mean flow and turbulence transport 
equations results in the following forms of the flow 
field governing equations: 


Continuity equation: 


3 V. 


aF 0 t 

2 '-5T*TT* f -° 


»„ere V, . if ♦ ( 2t ) 0 ' 5 f 


(ID 


Momentum equation: 


2SF 


ll + v It = 2 t(l - F 2 ) 3U e + a"F + 




t 3n 


3£ 


an 


2 )°.5 |S 
' 3n 


(12) 


Enthalpy equation: 

25F TT +v tS = -Tir + ^t7 +(25) 


3T 


0.5 

3n 


(13) 


Turbulence equations: 


aR 3R 4R F 3U U a / n , 3R 

2r F -i-+V — — = -2 e — r S. + JL i_ aR 0 ' 5 — X 

?r 3C t 3n 4 U 3C V 3n \* 3n 


3 2 R 0 5 sF (R - 3 R x) vR 0 - 5 

+ —J- + 2S(2c) U,b if + 2s , — — 

z 3 n o au 


2 ^ R x - 2 


aF 

2v 

a£ 

an 2v 2 \ 

9F 

" U 

3n 

3X ■ ,,2. 2 


(14) 


aR. 


2? 


3R. 


vR 0.5 (R - 3R y ) 
AU 3 


2tR. 


r 3R 1 

i 3 2 R 

5 y 1 
3n j 

lv 

2 if 3 

. 2v_3f J 


U e 3n " 


u > 
e i 


( 15 ) 


3 R 3R 7 2R F all U , / n , aR , 

2^F — ^ + V — £ = -2r — I £ + _e a_ / R 0.5 

35 v t 3n ^ U 3 E v 3n l 3n 


3 2r , vR 0.5 (R - 3R ) 2 2 

Z + 2 ^ 3 " 2 ^ R z TXT <“> 


3n 


U e 


2 rF + V ll _ _ 2t 2 FS + 3 _ 3 _ / R 0.5 *s\ 

? 35 t 3n " U e 3? 3 v 3n | AR 


+ f4 + (25)°- 5 R y f + 2 cS - 


vR°’ 5 2v 2 \ ,,,> 

U e A ~77J (1?) 


and 


39 x „ 39 o. Fe 3l e ^ , ”e 3 /. n 0.5 se 


2 CF + v = _2c L? '■ -S. + 3 -2-4- AR 1 

3 £ t 3n I 3£ v3r) 


3n 


+ — 0 + (20 0 * 5 R w + 2f;e ( + 


3n 


y 3n 


3C U 3n 3 x 
e 


vR 0 - 5 2v 2 ' 

U e A ” U 2 a 2 , 
e } 


( 18 ) 


where R = R x + R y + R z . 

The length scale a is required in the formulation 
of the above turbulence equations. This length scale 
distribution within the boundary layer is also assumed 
similar to the mixing length profile of a flat plate 
turbulent boundary layer. Within a small distance from 
the surface, the length scale is assumed to vary lin- 
early as function of the normal distance from the sur- 
face. Away from the surface, the length scale remains 
constant. This length scale distribution is shown in 
Fig. 2. The pivot point, y p , is equal to 10 percent 
of the theoretical laminar stagnation point boundary 
layer thickness 6. At the boundary layer edge, the 
length scale A e is related to the freestream 
turbulence intensity by A e = (0.25 - 3Tu)fi. 

To solve the above governing equations for the 
mean flow and turbulence properties, the initial and 
boundary conditions are derived analytically as 
fol lows: 

Initial conditions . The initial profiles of the 
mean flow and the turbulence correlations within the 
boundary layer at cp = cp 0 are determined from 
Eqs. (11) to (18). The following assumptions are 
made: 

1. At n 2> 


3R x 

1^=°* 


3R 

3R, 

% 

II 

o 

z 

3n 

if-o. 

and 


= o. 


I s = o. 

3n 
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In addition, the conditions of 


and 


F = 1, H = 1, S = c 5 (R x R ) 0 ' 5 , and e = 0 


x y' 


are also imposed at n 2 * An empirical constant, 

C 5 = 0.001 is used in this study. 

__ 2. The Reynolds normal stress components u 2 and 

v 2 are assumed to be independent of <p in the vicin- 
ity of the stagnation point. Therefore, at cp 0 , 

8 R 


3? 


= R 


(v*) 


(19) 


and 



3R v / -2 \ 

H y WZ) 

\ e o/ 

( 20 ) 

Similarly, 

aR 

TT = R z¥V 

( 21 ) 


ft= S W 

( 22 ) 

and 

f! = 

(23) 

are also assumed at <p 0 . 


3. The mean flow is symmetric about qp = 0 plane. 
In order to facilitate the analysis, these mean flow 
conditions are also assumed to be true at very small 
distances away from the stagnation point. Therefore, 
aF / a ^ = 0 and aH/as = 0 are imposed along the n 
direction at <? Q location. 

4. The conditions F = 0 , H = T w /T e , R x = 0, 

Ry = 0» R 7 - 0, S = 0, and e = 0 are also imposed 
at the cylinder surface. 

Applying the above assumptions to the turbulence 


transport equations, the following relations are found 
at Uo>n 2 ) location. 


D l. 5 
vR^ 
o,e 

HJ T 

o,e e 


x,o,e 


-2 


Re? 


4(dU 16 k ] 


U 


o,e 


o,e 


vFf" 

U ’a + 
o,e e 


2v 2 

ITT 

o,e e 


(24) 


vR 




1.5 

^e_ 


o,e A e 


y,o,e 


-2 


Recp" 


vR 


,0.5 

o,e 


U A 
o,e e 


2v 2 

U 2 A 2 
o,e e 


(25) 


f ! ( V 55 - 


2(du;dQ 0-e ^ 
’’o.e 


(26) 




2(dU/dc) 


o,e 


vR 0 ' 5 

o.e 

TJ T~ 

o.e e 


2v 2 

~2 7 

U o,e A e 


(27) 


where R Q p = 2 vaK^/ U§ e . Since 0 = 0 is imposed 
at n 2 * a§/dc = 0 occurs at the boundary layer edge. 

This zero gradient condition is also used for the 
boundary layer calculation at <p 0 . It follows from 
the assumption, Eq. (23), that 

f 3 U 0 ) - 0 (28) 

With the aid of fi(<p 0 ), f2(<p 0 )» and f 3 ( < Pp)» 
assumption 2 defines the mean flow streamwise turbu- 
lence gradient terms. Substituting Eqs. (19) to (23) 
into Eqs. (11) to (18) results in a set of ordinary 
differential equations with the boundary conditions 
specified at the surface (n = 0) and at n2* These 
equations are solved numerically for^the boundary layer 
flow property profiles at <p 0 = 0.04° in order to 
avoid the singularity at the stagnation point, <p = 0 . 0 °. 
These property profiles are the initial conditions for 
the downstream boundary layer flow field analysis. 

Boundary conditions . The turbulence double correl- 
ations along the boundary layer edge are also required 
to analyze the flow field. Turbulence kinetic energy 
and dissipation rate gradients in the normal direction 
are assumed to be zero (ak/ay = 0 , ae/ay = 0 ) at the 
boundary layer edge. The turbulence kinetic energy 
along the streamwise direction is determined by using 
the k-e two-equation turbulence modeling. The tur- 
bulence kinetic energy is then substituted into the 
Reynolds stress equations to determine the turbulence 
double correlations along the boundary layer edge. 

Introducing Eqs. (1) and (4) into the k-e equa- 
tions and seeking the following forms of solutions: 

k = vaK(c) and e = va 2 E(c) (29) 

the k and e turbulence modeling equations become 



From Eqs. (5) to (10) the initial conditions: 


K - K 2 and E = E 2 (32) 


are defined at ( s 2 , n 2 ) • 


With the assumption that the gradients aR x /an, 
aRy/an, and aR z /9n vanish at the boundary layer 
edge, Eqs. (16) to (18) give the following relations 
among the mean velocity, the Reynolds normal stresses, 
and their mean flow streamwise gradients; 


dR x,e 

d£ 


4R dU 


U. 


d£ 


rT 5 vR° -5 R v . , 2 

v e e x,e 2v n 

IL“Ta 7~ "X? x,e 


e e 


U a*“ 
e e 


(33) 
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“Ve ” R e' 5 » "“''Ve 2. 2 . 

-77V, 

e e 


(34) 


and, 

dR 

~W 


z.e 


2R 


dll 


vR 1 - 5 vR 0 - 5 


2v R 

uR z,e 


(35) 

The above equations are solved for the Reynolds 
normal stresses. The sum of the Reynolds normal 
stresses R e should be consistent with the turbulence 
kinetic energy obtained from Eqs. (30) and (31) (see 
Appendix A) . 

The previous analysis of the initial conditions 
for Region 3 defines the Reynolds normal stresses at 
(^2^2)* These initial conditions are used to inte- 
grate Eqs. (33) to (35) along the downstream direction. 
The Fourth Order Runge-Kutta numerical method is used 
for the integration. The turbulence kinetic energy 
obtained previously from Eqs. (30) and (31) is used to 
evaluate the derivatives. The integration step is 
adjusted until the turbulence kinetic energy 

(R x e + R y e + R z e) u e^ satisfactorily converges 
to the corhesponalng value predicted by the k and e 

equations. The corresponding Reynolds normal stresses 
are the boundary conditions at the boundary layer 
edge. 

Finally, the following conditions are also 
imposed: 


1. at the surface n = 0, 

F = 0, H = T w /T e , R x - 0, R y = 0, 

R z = 0, S = 0, and e = 0 

2. at the boundary layer edge n = n2> 

F e = 1, H e =1, Sg = C5( Rx , e^y ,e)^‘ 

and e=0 

A numerical computational procedure to calculate 
the boundary layer flow field described by Eqs. (11) 
to (18) was developed in a previous study (12). This 
existing numerical computational procedure is used to 
calculate the turbulence correlations and the mean 
flow properties in Flow Region 3. 

RESULTS AND COMPARISON WITH EXPERIMENT 

In order to verify the above theoretical analysis, 
a numerical computation was performed to compare the 
analysis with existing experimental measurements. A 
set of experimental data on the effect of freestream 
turbulence on heat transfer from heated cylinders 
placed normal to airstream was reported in Ref. 1. 

The two point turbulence longitudinal correlations in 
the freestream were also measured in that experiment. 
These data are used to determine the turbulence longi- 
tudinal microlength scale which is required to start 
the present numerical calculations. Thus, the test 
conditions of this reference were used as input to the 
present numerical computational procedures. The tur- 
bulence longitudinal microlength scales are obtained 
by curve-fitting the experimental data using Eq. (3). 
Based on the freestream mean velocity, the freestream 
turbulence intensity, and the turbulence longitudinal 


microlength scale, the turbulence kinetic energy, the 
turbulence double correlations, and the mean flow pro- 
perties within the stagnation flow field are computed. 
The surface heat transfer rate is also calculated from 
the mean temperature distribution. In the following 
sections these numerical results are described and are 
compared with the experiment of Ref. 1. 

Turbulence 

To show the variation of the turbulence kinetic 
energy along the stagnation point streamline, some 
analytical results are plotted in Fig. 3. Except for 
the lowest turbulence case, the turbulence kinetic 
energy increases within the inviscid flow region in 
the direction toward the surface and reaches a maximum 
value at the boundary layer edge. The turbulence 
kinetic energy then decreases continuously and van- 
ishes at the cylinder surface. This turbulence 
kinetic energy variation is similiar to the results 
in Ref. 10. However, the present analysis also pre- 
dicts a reduction in the turbulence kinetic energy 
along the stagnation point streamline when the free- 
stream turbulence intensity and Reynolds number are 
small, (Tu = 0.013 and Re = 1.90xl0 5 ). With con- 
stant Reynolds number, increasing the freestream tur- 
bulence intensity also increases the turbulence 
kinetic energy at the boundary layer edge (n ^ 5). 

Examples of the Reynolds normal stress profiles 
within the boundary layer at the stagnation point are 
shown in Fig. 4. These results are for the cases with 
Re = 2.5xl0 5 and Tu = 0.017 and 0.028. The Reynolds 
normal stresses decrease continuously from the bound- 
ary layer edge in a direction normal to the surface. 

The Ry value was predicted to be the largest for a 
small region very near the surface. These Reynolds 
normal stress profiles differ significantly from those 
of a flat plate turbulent boundary layer (18) where 
the largest Reynolds normal stress was foun8" in the 
x-direction. 

Analytical predictions of the mean velocity com- 
ponent, F = U/U e within the boundary layer near 
the stagnation point are shown in Fig. 5. These pro- 
files are compared with the mean velocity profile of 
the theoretical laminar boundary layer. The compari- 
son shows that increasing the freestream turbulence 
increases the velocity near the surface and decreases 
the velocity near the boundary layer edge. Existing 
measurements ( 19 ) of this velocity component F with 
Re = 2.5xlCr and Tu « 0.05 are also plotted in this 
figure. The analysis agrees with the measurements. 

Example of the computational results of the tur- 
bulence kinetic energy and the Reynolds normal 
stresses along the boundary layer edge are shown in 
Fig. 6. In this figure, the turbulence quantities 
are nondimensionalized with the corresponding values 
at (S2»n2)* The turbulence kinetic energy increases 
continuously along the downstream direction. Large 
increase in the turbulence kinetic energy is found for 
the case with small freestream Reynolds number and 
turbulence_i ntensi tv. The Reynolds normal stress com- 
ponents, Ue and Ve, increase monotonically at 
the downstream locations. However, the wjT compon- 
ent shows a reduction near the stagnation point. The 
effect of the freestream turbulence longitudinal micro- 
length scale on the Reynolds normal stresses along the 
boundary layer edge is shown in Fig. 7. Increasing 
the turbulence longi tudi naljuicro length scale is shown 
to decrease the value of u e^ u o,e and increases 
the value of w|/Wq e along the boundary layer 
edge. No significant difference in the value of 
vf/v^e is calculated. 
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Surface Heat Transfer 

The present analytical results of the surface heat 
transfer is defined in terms of the Nusselt number; 

" u ■ (w) y ,„/ T » - V <36 > 

The mean temperature gradient (3T/ay)y = o at the sur- 
face is calculated numerically from tne mean temper- 
ature profile, which in turn, was obtained previously 
from the solution of the boundary layer flow. 

The analytical results for the stagnation point 
surface heat transfer rate over a range of the free- 
stream Reynolds number and the turbulence intensity 
are compared with experimental data U) and existing 
correlations U,2) in Fig. 8. The surface heat trans- 
fer rate is expressed in the form of the Froessling 
number (NuxRe“^* 5 ). Linear variation in the turbulence 
longitudinal microlength scale as function of the 
freestream Reynolds number is also assumed in the 
present calculation. With constant freestream turbu- 
lence intensity, the longitudinal turbulence micro- 
length scales within ±10 percent of the values at 
Re = 2.5x10^ were used in the numerical computation 
for a range of 1. 5£Rexl0~ 5 <3.0. Including the effect 
of the freestream turbulence microlength scale, the 
present analysis calculates a larger stagnation point 
surface heat transfer rate than is indicated by the 
correlations in Refs. 1 and 2. 

Examples of the computational results of the tur- 
bulence double correlation -vt/U e Te within the 
boundary layer at the stagnation point are shown in 
Fig. 9. The maximum turbulence correlation occurs 
near the surface. Increasing the freestream turbu- 
lence Jjitensity also increases the turbulence correla- 
tion vt. Together with the characteristics of the 
Reynolds normal stress component Ry as shown in 
Fig. 4, Eqs. (13) and (18) predict a high heat trans- 
fer rate at the surface. It is postulated that the 
freestream turbulence penetrates into the near sur- 
face flow region and induces large turbulence cor- 
relation vt. Thus, a surface heat transfer rate 
larger than that of a laminar flow occurs. 

Some analytical results of the surface heat trans- 
fer rate within a small distance from the stagnation 
point are plotted in Fig. 10 in terms of Froessling 
number. Approximately constant Froessling number is 
calculated only for the case with Tu = 0.013. For 
Tu = 0.017 and 0.028, the analytical results show a^ 
reduction in the Froessling number near the stagnation 
point. With high freestream turbulence intensity, a 
large turbulence dissipation rate may occur near the 
stagnation point (Eq. (2)). This turbulence dissi- 
pation rate reduces the rate of the increase in the 
turbulence kinetic energy. This rate of increase in 
Ke may not be sufficient to induce a constant sur- 
face heat transfer rate. 

To investigate the effect of the turbulence longi- 
tudinal microlength scale on the surface heat transfer 
rate, the analytical results of the turbulence kinetic 
energy and the Froessling number for the cases of 
X p / D = 0.0053 and 0.0079 with Re = 2.5x10^ and 
Tu = 0.028 are compared in Fig. 11. With 
xp/0 = 0.0079, the present analysis computes the 
experimental data in Ref. 1. A lower level of the 
turbulence kinetic energy and a smaller Froessling 
number are also predicted for the case with small 
freestream turbulence longitudinal microlength scale. 
These computational results also show a reduction in 
the Froessling number near the stagnation point. 

In this study, the approximation, 

S e = C 5 (R X s e Ry >e )^* , is used along the boundary layer 


edge. This approximation is similar to the the cor- 
relation observed from the flat plate turbulent bound- 
ary layer study (_18) . The constant c5 = 0.001 is 
found to give the best computational results in the 
present analysis. In an earlier study (12) , a dif- 
ferent value of C 5 = 0.0001 was used because the 
boundary layer edge location and the corresponding 
turbulence kinetic energy (or the Reynolds normal 
stresses) were not analyzed theoretically. If one 
uses the small value of cs, then a small (nonzero) 
uv is imposed along the boundary layer edge. This 
approximation is different from the k-e two-equation 
turbulence modeling, which requires uv to be zero 
for the external inviscid flow field in this study. 
However, the requirement of uv = 0. is itself an 
assumption. 

The assumption -vt =0 is also imposed at the 
boundary layer edge. Intuitively, the turbulence 
transport equations can also be used to study the_ 
variation of the turbulence double correlation _vt 
within the external inviscid flow region. If vt i_s_ 
nonzero in the external flow field, then different vt 
may occur within the boundary layer flow region and 
the present analysis will predict different stagnation 
flow surface heat transfer rate. If this is true, it 
might resolve the different level of the stagnation 
flow surface heat transfer rate reported in various 
experiments. Although, the task requires the vt 
boundary condition within the freestream, the effect 
of the freestream temperature fluctuation on the sur- 
face heat transfer rate of the stagnation flow can be 
studied using the present analysis. 

CONCLUSIONS 

Theoretical analysis was formulated to model the 
effect of the freestream turbulence on the surface 
heat transfer of a stagnation flow region. The flow 
field of interest is around the forward stagnation 
point of a circular cylinder in a turbulent flow. The 
k-e two-equation turbulence modeling was used to 
analyze the turbulence kinetic energy within external 
inviscid flow. The time-averaged turbulence transport 
equations were employed to predict the turbulence 
double correlations within the boundary layer. Based 
on the boundary layer edge conditions of the turbu- 
lence kinetic energy, the transport equations of the 
mean flow and turbulence are solved numerically for 
the turbulence double correlations, the mean velocity 
and the mean temperature. 

The present analytical results show that, depend- 
ing on the freestream turbulence properties, different 
levels of turbulence kinetic energy occur along the 
boundary layer edge. These different levels of turbu- 
lence kinetic energy induce different surface heat 
transfer rates. Increasing the freestream turbulence 
intensity or the turbulence longitudinal microlength 
scale also increases the stagnation region surface 
heat transfer rate. It is postulated that the free- 
stream turbulence penetrates into the boundary layer 
and induces high surface heat transfer rate. 

The analysis presents a procedure to calculate the 
stagnation point surface heat transfer rate as func- 
tions of the freestream turbulence intensity, turbu- 
lence longitudinal microlength scale and Reynolds 
number. Due to limited data on the turbulence longi- 
tudinal microlength scale, this analysis was verified 
only for some cases with Reynolds number of the order 
of 1(P and the freestream turbulence intensity of less 
than 5 percent. Considering the freestream turbulence 
longitudinal microlength scale, this analysis predicts 
a higher stagnation point surface heat transfer rate 
than was obtained from the existing correlations 
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relating the freestream turbulence intensity, Reynolds 
number and the surface heat transfer rate. 

APPENDIX A 

The turbulence kinetic energy and its dissipation 
rate along the stagnation point streamline and the 
boundary layer edge are calculated respectively from 
Eqs. ( 6 ) and (7), and Eqs. (30) and (31). Finite dif- 
ference method is used to solve these equations for 
the K and E distributions which satisfy the bound- 
ary conditions. 

The dimensionless Eqs. ( 6 ) and (7) can be written 
as 

A nl = A n 2 6 + A n 3 ^ ^ + A n 5 (A1 ] 

with G = K or E 

Similarly, Eqs. (30) and (31) can also be written 
as 

\lf = \2 G + \3^(\4f) + \5 < A2 > 

with G = K or E. 

In the above equations, A^-j and A n i may be con- 
stant, or the functions of the turbulence properties, 
or coordinates s and n. 

A similar finite difference scheme is used to 
convert Eqs. (Al) and (A2) into finite difference 
forms. The derivatives at each grid point m are 
approximated by central finite difference. A finite 
difference equation 

B l^m + l + B 2 G m + B 3 G m-l = B 4 ( A3 ) 

is found at each grid point m. The coefficients 
B-j (i = 1 to 4) are functions of the various coef- 
ficients or A^-j , and the grid size at the grid 

points m+ 1 , m, m- 1 . 

Applying the previous mathematical manipulation 
at every grid point within the interval nl > n _> n 2 
or $2 _> Z 1 £3 results in a set of finite difference 
equations for Eq. (Al) or Eq. (A2). The boundary con- 
ditions of G are specified at nl and n2 ( or £2 
and £ 3 ). This set of equations are solved for the G 
values at each grid point m. The grid size within 
the computational interval is generated using the model 
in Ref. 20. Two hundred to four hundred grid points 
are used in the numerical computation. 

Equations ( 8 ) to (10) describe the boundary con- 
ditions at n = nl and n2» However, the value 
of n2 not yet defined. Therefore, an arbitrary 

value is assigned to n 2 to start the numerical cal- 
culation. The existing analysis (_10) indicates that 
n 2 is approximately equal to 5.0. Linear variations 
in K and E are first assumed. The coefficients of 
the finite difference equations are evaluated using 
these K and E distributions. Thus, two sets of 
linear algebraic equations are obtained to relate K 
and E at the grid points. Since K and E are 
prescribed at nl and n 2 > these equations are solved 
to update the K and E distributions using the 
method of successive substitution. The above proce- 
dures are repeated until K and E converge at each 
grid point. The boundary values of K and E at 
n2 are adjusted to assure that K and E profiles 
satisfy the boundary conditions at nl and n 2 * 

The convergence in K and E is sensitive to 
their assumed values at the start of the computations. 
Based on the k-e two-equation turbulence modeling. 


Strahle (11) developed analytical links between the 
turbulence properties k and e at the stagnation 
point and in the freestream. Neglecting the diffusion 
terms and using the same value for the empirical con- 
stants ci and C2, (ci = C 2 = 1.90), these analyti- 
cal links can be written in nondimensionalized forms 
as 


k 


2 “ 



1.111 


k 


2.111 

F 


(A4) 


and 


E r (0.6 K p ) 


2 .111^-1/1. Ill 


(A5 ) 


With the isotropic assumptions, Eq. (2), the 
approximate boundary conditions for K and E at n2 
are calculated from the previous relations. These K 
and E values are used to start the numerical 
computations. 

The previously mentioned numerical calculation 
provides the conditions of K and E at £ 2 * Based 
on these conditions, the K and E along the bound- 
ary layer edge are also calculated from the sets of 
finite difference equations corresponding to Eq. (A2) . 
The computational domain, S2 1 £ <_ £3 limited to a 
small region near the stagnation point, 

(0.04° < cp < 15.0°). 

First, K and E at £3 location are estimated 
and linear distributions in K and E between ££ 
and £3 are assumed. The coefficients of the finite 
difference equations are calculated using the assumed 
K and E distributions. Two sets of linear algebraic 
equations relating the K and E at the grid points 
are obtained. These algebraic equations are solved 
with the method of successive substitution to update 
the K and E distributions. This iterative com- 
putational procedure is repeated until K and E 
converge at each grid point. The values of K and 
E at S 3 are also adjusted in each iteration to 
facilitate the convergence in the numerical 
computations. 
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DISTANCE NORMAL TO SURFACE, y/6 



FIGURE 2. - MODELING LENGTH SCALE OF THE TURBULENCE 
CLOSURE ASSUMPTIONS. 
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FIGURE 3. - TURBULENCE KINETIC ENERGY ALONG THE STAGNATION 
STREAMLINE. 




FIGURE A. - REYNOLDS NORMAL STRESS PROFILES WITHIN THE 
STAGNATION POINT BOUNDARY LAYER, Re = 2.5x10 5 . 
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FIGURE 5. - MEAN FLOW STREAMWISE VELOCITY COMPONENT 
WITHIN THE STAGNATION POINT BOUNDARY LAYER. 

Re = 2.5x10 s . 
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FIGURE 6. - BOUNDARY LAYER EDGE TURBULENCE KINETIC ENERGY AND CORRESPONDING FROESSLING NUMBER, Tu = 0.013. 




FIGURE 7. - EFFECT OF THE TURBULENCE LONGITUDINAL MICROLENGTH SCALE ON THE REYNOLDS NORMAL 
STRESSES ALONG THE BOUNDARY LAYER EDGE, Re = 2.5x10 5 AND Tu = 0.028. 
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STAGNATION-POINT FROESSLING NUMBER, Nu Re 



FIGURE 8. - COMPARISON OF THE STAGNATION POINT 
FROESSLING NUMBER BETWEEN THE ANALYSIS AND EXPERI- 
MENTAL CORRELATIONS. 
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TURBULENCE CORRELATION, -vt/U e T e 


FIGURE 9. - PROFILES OF THE TURBULENCE CORRELATION vt 
WITHIN THE STAGNATION POINT BOUNDARY LAYER. 
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FIGURE 10. - TURBULENCE KINETIC ENERGY AND REYNOLDS NORMAL STRESSES ALONG THE BOUNDARY LAYER EDGE. 
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FIGURE 11. - EFFECT OF THE TURBULENCE LONGITUDINAL MICROLENGTH SCALE ON THE TURBULENCE KINETIC ENERGY AT THE 
BOUNDARY LAYER EDGE AND FROESSLING NUMBER, Re = 2.5x10 5 AND Tu = 0,028. 
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